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SUMMARY 

A nonlinear laminate analysis is described for prc.iicting the mechani- 
cal behavior (stress-strain relationships) of angleplied laminates in which 
the matrix is strained nonlinearly by both the residual stress and the mech- 
anical load and in which additional nonlinearities are induced due to pro- 
gressive fiber fractures and ply relative rotations. The nonlinear laminate 
analysis (NLA) is based on linear composite mechanics and a piece-wise 
linear laminate analysis to handle the nonlinear responses. Results ob- 
tained by using this nonlinear analysis on boron-fiber/aluminum-matrix 
angleplied laminates agree well with experimental data. The results shown 
illustrate the in situ ply stress-strain behavior and synergistic strength 
enhancement. 


INTRODUCTION 

Metal matrix composites, especially boron-fiber/aluminum matrix, have 
been investigated intensively in recent years and are now considered seri- 
ously for use in aerospace structures (Refs. 1 and 2). Predicting the mech- 
anical behavior and structural response of components fabricated from metal 
matrix composites requires use of various mathematical models. These 
models, for example, relate stresses to applied forces, stress intensities 
at the tips of cracks to nominal stresses, buckling resistance to applied 
force, and/or vibration response to excitation forces. The Models just 
mentioned require initial tangent and strain-dependent stress-strain rela- 
tionships. Experimental data indicate that the stress-strain curves of uni- 
directional composites are: (1) slightly nonlinear along the fiber direc- 

tion near fracture; (2) mildly nonlinear transverse to fiber direction; and 
(3) extensively nonlinear in intralaminar (in-plane) shear. In addition, 
analysis of experimental data shows that the lamination residual stresses 
(the fabrication process induces thermal strains) may be of sufficiently 
high magnitude to strain the aluminum matrix nonlinearly in certain ply 
orientation configurations (Refs. 3 to 5). When the matrix is strained 
nonlinearly, the stress-strain relationships of the laminate may become 
load-path dependent. Though the nonlinear response of composite behavior 
has received some attention (Refs. 5 to 7), the lamination residual strain 
nonlinearity, including its effects on subsequent loading, has not been 
examined. It is the purpose of tnis investigation to describe a nonlinear 
laminate analysis for predicting the nonlinear mechanical behavior (stress- 
strain relationships) of angleplied laminates in which the matrix is 
strained nonlinearly by both tne residual stress and the mechanical load. 
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The nonlinear laminate analysis (NLA) is an approximate analysis. It 
consists of linear composite mechanics and a piece-wise linear laminate 
analysis to handle the nonlinear response of the in situ matrix and individ- 
ual fiber fractures. The NLA, as used in the investigated described herein, 
has no programmed provisions for the automatic detection of laminate frac- 
ture. This provision will be incorporated in the near future. 

THEORETICAL BASIS 

The computational procedure developed for the nonlinear laminate analy- 
sis (NLA) of boron/ aluminum angleplied laminates is approximate. It con- 
sists of linear composite micromechanics and piece-wise linear laminate 
analysis to handle the cumulative effects of the nonlinear response of the 
in situ matrix. Briefly, tangent properties of the current cumulative 
matrix strain level are used to compute ply properties using composite 
micromechanics. These ply properties are then used to determine laminate 
properties using linear laminate theory. The load is incremented and the 
ply stresses/strains are calculated using linear laminate analysis and the 
previous increment laminate and ply properties. The ply stresses and 
strains are used to calculate matrix strains in conjunction with ply and 
matrix previous increment properties. The matrix strains thus calculated 
are amplified using integrated matrix magnification factors (Ref. 5). The 
amplified strains are added to the cumulative strains and the matrix tangent 
properties are determined using "table look-up" for the next load incre- 
ment. At this point also the percent of fibers fractured, if any, is deter- 
mined and is used to modify the longitudinal ply properties. 

The residual stresses due to a temperature difference between process- 
ing and use temperature are determined in slightly different way. The 
macro- and micro-residual strains are calculated using linear composite 
mechanics and average matrix thermal and mechanical properties. These aver- 
age properties are used as inputs at the micromechanics analysis level to 
obtain the ply properties. The ply properties and the temperature differ- 
ence between processing and use-tempe>^atures are used as inputs as the lami- 
nate analysis level to compute the residual strains and stresses in the 
plies. The integrated residual strains in the matrix are computed using the 
ply strains and the strain magnification factors (Ref. 5). This approach 
predicted results which were in good agreement with measured data for ini- 
tial tangent properties (Ref. 5). 

Two key assumptions form the basis for the development of this NLA. 
These are: (1) the in situ matrix is nonlinearly orthoropic; and (2) the 

orthotropic properties associated with each strain component are independent 
and can be determined from table look-up of normal and shear stress-strain 
data. These assumptions introduce some approximations in the stress-strain 
relationships (constitutive relationships) of the in situ matrix. However, 
these approximations are thought to be minimal and, therefore, to have 
negligible effect on the results because of: (1) the in situ matrix is 

highly restrained by the proximity of the fibers and no flow rule in the 
"classical plasticity sense" is operative; (2) the relatively large magni- 
tude of microresidual stresses in the matrix and (3) and most important, 
force equilibrium and strain compatibility are automatically satisfied 
individually and collectively at both the micromechanics and macromechanics 
levels. 

The computational procedure for NLA is by necessity in the form of a 
computer program. In the present investigation, the computer program was 
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generated by modifying an available computer program (Ref. 8). This program 
consists of a collection of functional modules (routines) which carry out 
composite micromechanics, macromechanics, and laminate theory. The inputs 
to this computer program are constituent material properties, laminate con- 
figuration, and composite geometry. 

A qualitative flow chart of the resulting computer program is shown in 
Fig. 1. The temperature difference is an input. Note the provision to 
check the positive definite condition of the array of the stress-strain 
relations of the aluminum matrix. Violation of this condition leads to pre- 
dictions of negative moduli of elasticity for both the ply and the compo- 
site. The asterisk indicates where the initial tangent properties for the 
laminate are oetermined. 

It can be seen in the flow chart that the NLA procedure is “tangent- 
marching." That is, the laminate response at the next load increment is 
predicting using previous increment tangent properties without local itera- 
tion (predictor but not corrector). Local iteration added only to the in- 
creased accuracy of the results but not to the conceptual formal structure 
of the NLA procedure. Besides, local iteration (predictor corrector ap- 
proach) can readily be incorporated as needed. 

The NLA used in the present investigation has several important fea- 
tures when compared with other approaches. These are: (1) composite micro- 

mechanics is used to describe the matrix nonlinear behavior; (2) progressive 
fiber fractures are tracked at the micromechanics level; (3) initial or 
residual strains are determined prior to mechanical load application; 

(4) development of suitable expressions for the approximate representation 
of the interfiber matrix strain variation; and (5) generally orthotropic 
nonlinear behavior of the in situ matrix is permitted for both initial 
strains and subsequent mechanical loadings. The main advantages of this 
approach are: (1) computer running time economy, (2) ease of input data 

preparation, (3) measured nonlinear stress-strain data for the ply under 
cc.T.oined load is not required, and (4) the three-dimensional tangent proper- 
ties for both ply and composite are computed routinely as a part of the 
analysis. 

As was mentioned in the introduction this NLA does not have provisions 
for automatically detecting fracture of the laminate. This is determined by 
the user from the ply stress or strain states and/or from the execessively 
large strains of the laminate. 


INPUT/OUTPUT 

The input data required for the nonlinear computational procedure con- 
sists of tabular data of matrix thermal properties, stress-strain data due 
to mechanical load, and percent of fibers fractured due to mechanical load. 
The thermal data of the matrix was obtained from Ref. 9 and are shown in 
Table 1. 

The mechanical load stress-strain data were obtained from axial tension 
and from torsion tests on thin tubes made from 6061-0 aluminum.. Briefly, 
the tube test section was 2 inches inside diameter, 10 inches long, and 
0.060 inch in wall thickness. The tube was machined to the test section 
dimensions from commercial tube stock allowing 1 inch at each end for 
gripping. The tubes were instrumented at their midlength with strain gage 
rosettes. The tubes were loaded in a multiaxial testing machine under ten- 
sion and torsion loads. Load and strain data were recorded at regular load 
intervals. These data were reduced to engineering stress-strain data using 


3 


a strain-gage data-reduction program (Ref. 10). Tensile load data suitable 
for input into the computational procedure are shown in Table 2. The corre- 
sponding shear data are shown in Table 3. The data for percent of fibers 
fractured was obtained from Ref. 11 and is summarized in Table 4. 

Additional input data required include boron fiber material properties, 
composite geometry laminate configuration, load and load increment, tempera- 
ture difference, and the correlation coefficients described in Ref. 8. The 
fiber properties used were: modulus = 57 x 1Q6 psi, Poisson's ratio = 0.2, 

and thermal coefficient of expansion = 2.8x10“° in/in"F. This input data 
set is general for boron/ aluminum angleplied laminates with boron fiber in 
6061-0 aluminum alloy. 

A large amount of information is generated during the solution pro- 
cesses and printed out in addition to the input data. This information con- 
sists of the following general sets: (1) All the input data is printed out 
in table form for ease of cross checking and reference; (2) The results of 
the lamination residual stresses are printed out. These results include ply 
stresses and strains, matrix stresses tid strains, and matrix moduli and 
Poisson's ratios corresponding to t^, ."esidual stress/strain; (3) First load 
increment results are printed for initial tangent moduli and Poisson's 
ratios for both the ply and the laminate (Ref. 5); (4) Any or each load 
increment results are printed at user request. These results include ply 
and laminate tangent properties for that load increment and cur.ulative 
properties up to that load increment. All these results are printed out 
automatically at the last load increment; (5) Final cumulative results are 
printed after the last load increment. These are in table form listed at 
each load increment and include ply stresses and strains, ply equivalent 
matrix strain, ply relative rotation, laminate stresses and strains and 
laminate geometry changes for length, width and thickness. The final cumu- 
lative results can be used to plot nonlinear stress-strain curves for the 
laminate, in situ ply stress-strain curves, ply stress/strain influence co- 
efficients and interply relative rotation. Some of these will be illus- 
trated later in the discussion of results. 

RESULTS, COMPARISONS AND DISCUSSION 

The nonlinear laminate analysis (NLA) described previously was used to 
predict the stress-strain behavior of boron/aluminum (B/A) angleplied lami- 
nate (APL) with generic laminate configuration [02/*©]$ loaded in ten- 
sion at various angles to the 0* ply direction. The laminates were made 
from 4-mil diameter boron fiber in a 6061-0 aluminum alloy. The fiber vol- 
ume ratio was about 0.50. The laminate fabrication and experimental program 
to generate the data are described in detail (Refs. 3 to 5). 

The temperature difference used in the residual strain computations in 
this investigation was 700* F. The reason for this choice is tha| the 
matrix starts supporting stress greater than 1000 psi between 800 and 
700* F (Table 1). Also, this temperature difference was used to predict the 
initial tangent properties (Ref. 5), which correlated with measured data. 

The load at fracture was applied to the laminate in twenty equal incre- 
ments. The computer running time was about 4 seconds for the twenty incre- 
ments in the UN I VAC 1110. 

NLA results are comgared^with measured data in Fig. 2, for a unidirec- 
tional laminate loaded 0 , 10*, and 90* to the fiber direction. The results 
plotted are laminate stress versus laminate strain. The agreement is ex- 
cellent except for the final strain in the 10* off-axis specimen. Part of 
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this difficulty may be due to the extrapolation method that is used to 
determine the measured strain at fracture (Ref. ?.0). The results in Fig. 2 
show that unidirectional boron/aluminum laminates exhibit: (1) mild non- 

linearity when loaded parallel (longitudinal) or perpendicular (transverse) 
to the fiber direction; (2) extensive nonlinearity when loaded at 10* off- 
axis, and (3) relatively low strain (0.08 percent) and stress (16 ksi) to 
fracture in the 90* (transverse) direction. This fracture strain is only 
about one-fifth of that of the aluminum alloy (5 percent) and the fracture 
stress is about 130 percent of the corresponding stress of the aluminum 
(Table 2). The longitudinal fracture stress of about 180 ksi coincides with 
approximately 7 percent- of fibers fractured (interpolated) Table 4. 

The 10* off-axis stress-strain curve is shown here because it is used 
to determine the intralaminar shear (in-plane) stress-strain curve 
(Ref. 12). The result is shown in Fig. 3. It is included here to indicate 
the large intralaminar shear strain to fracture of about 6 percent which is 
several times larger than that in Fig. 2. The intralaminar shear strain at 
fracture is about six times greater then the “yield" shear strain of the 
aluminum alloy (1 percent). Table 3. The intralaminar sheaf stress at frac- 
ture is about 12 ksi which is about 150 percent of the shear yield stress of 
the aluminum alloy. 

NLA stress-strain results for a [02/*15]5 APL loaded at -80* to the 
0* ply direction are compared with measured data in Fig. 4. Again, the 
agreement is very good except for the predicted strain at fracture which is 
about ten times the measured value. In this case it may be considered to 
imply that the laminate fractured, since a relatively large strain (about 
4 percent) was produced by the last load increment. 

The in situ 0* ply stress-strain behavior predicted by the NLA is 
plotted in Fig. 5. The schematics in this figure depict the loading in the 
laminate, the laminate configuration, the ply orientation, and the stresses 
plotted. The transverse stress ( 0122 ) curve shows extensive nonlinear 
behavior. This stress is about 20 ksi near fracture. The corresponding 
strain is about 5 percent. 

This 20 ksi stress is about 25 percent greater than the uniaxial value 
(16 ksi). The corresponding strain is about 50 times greater than the uni- 
axial value (0.1 percent). These comparions show significant in situ en- 
hancement or “synergistic effect" especially for the transverse strain. The 
longitudinal stress (om) is compressive and is about 40 ksi near frac- 
ture. The intralaminar shear stres| (o]i2) is only about 2.5 ksi. Com- 
parable in situ results for the +15* ply are shown in Fig. 6. The curves 
are similar for ^o]22 and o-jn with those of the 0* ply (Fig. 5). 

However, the +15* ply has an intralaminar shear stress of about 21 ksi and a 
corresponding strain of about 0.6 percent. The intralaminar shear stress is 
almost twice that of the uniaxial value (12 ksi) Fig. 3. Most of a-\i 2 
is due, to lamination residual stress since it remains practically constant 
with strain. The shear strain on the other hand is only about 10 percent of 
the uniaxial value. These results show that the [02/*l§]s B/A AgL 
failed by transverse tensile ply fracture in both the 0 and +15* plies as 
would be expected. However, the significant in situ enhancement synergistic 
effect indicated by the NLA results has not been reported previously even 
though some part may be expected. In addition the magnitude of this 
“synergistic effect" can only be predicted by NLA based at the composites 
micromechanics level. 

The 0* ply cumulative stresses are shown in Fig. 7 to illustrate the 
redistribution of stresses as one of the ply stresses becomes excessively 
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nonlinear. Note the longitudinal stress (om) starts increasing rapidly 
at a laminate stress of about 14 ksi as the transverse stress (0122) 
approaches a horizontal tangent and therefore cannot support additional 
stress. The intralaminar shear stress (o]i2) continues to remain 
negligible. 

NLA results are compared with measured data in Fig. 8 for the same 
[02/*15]5 B/AL APL Lut loaded at 5* to the ''0'' ply direction. See sche- 
matic. In this figure results are compared for axial stress (ocxx)» 
Poisson's strain (e^yy), and "coupling shear" strain (e^xy) versus axial 
strain (e^xx)* comparisons are in very good agreement. Again, there 
is some discrepancy for the last load increment strain. The reason for this 
discrepancy is the extrapolation method used to determine the measured 
strain at fracture^as was mentioned previously. 

The in situ 0* ply stress-strain behavior predicted by the NLA is 
plotted in Fig. 9. As expected, the longitudinal stress (o-qi) dominates 
and reaches about 160 ksi at fracture. This approaches the uniaxial longi- 
tudinal fracture stress of 182 ksi. The o^her two stresses are relatively 
small. Corresponding stresses for the +15 ply are plotted in Fig. 10. The 
intralaminar shear stress (0^12) starts at about -25 ksi (which is the 
value of lamination residual stress) and decreases to about 20 ksi. The 
transverse stress starts at about 10 ksi and decreases to about 5 ksi. The 
ply transverse strain reverses sign from initially positive (due to residual 
stress) to negative due to the applied stress. The in situ -15" ply stress- 
strain behavior is plotted in Fig. 11. The transverse stress (0^22) 
starts at about 8 ksi (residual stress) and decreases to about 3 ksi indi- 
cating that this stress is predominantly residual. The intralaminar shear 
stress (0112) starts at about 17 ksi and decreases to about 10 ksi. The 
corresponding strains reverse sign from positive to negative. The -15’ ply 
cumulative stresses are plotted in Fig. 12. The longitudinal stress 
(°lll) starts from slightly negative and remains linear to fracture with 
increasing laminate stress. The transverse (0122) and intralaminar shear 
(0112) stresses vary only slightly with laminate increasing stress. 

The NLA results for this laminate indicate that this laminate failed by 
longitudinal tension controlled mainly by fiber fractures. The other two 
stresses indicate negligible effect on the laminate strength. 

The comparisons described previously were selected to illustrate the 
versatility of the NLA. In addition they were selected to illustrate lami- 
nate material nonlinear behavior, the effects of lamination residual 
stresses and in situ synergistic effects. Taken collectively the results 
show that the NLA described provides an effective analysis capability to 
study and determine the material nonlinear behavior of boron/ aluminum APL. 
Such a NLA should also be applicable to metal matrix composites in general 
provided the required input data is available. The NLA has the capability 
to handle other membrane loads and bending loads and combinations. This 
part of the capability, however, has not as yet been validated with measured 
data. 

Improvements to this NLA capability include local convergence, dynamic 
load increment and graphical display of predicted results. Additional 
capability includes temperature incrementation for the residual stress, 
fracture detection criteria, and combinations of thermal and mechanical 
loads incrementation and the appropriate convergence schemes. 
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CONCLUSIONS 


A nonlinear laminate analysis (NLA) for predicting the mechanical be- 
havior of angleplied metal matrix laminates has been formulated and pro- 
grammed. The NLA was formulated using composite micromechanics, composite 
macromechanics, and laminate theory. This NLA allows the in situ matrix to 
be nonlinearily orthotropic and accounts for residual stresses, progressive 
fiber fractures, and relative ply rotations. Results predicted by this NLA 
are in excellent agreement with experimental data from boron-fiber/aluminum 
matrix composites except for the final strain in some cases. These results 
show that in situ synergistic effects enhance significantly the uniaxial 
composite transverse and intralaminar shear strengths. These results also 
show that a NLA provides an effective analysis capability to study and 
determine the material nonlinear behavior of metal-matrix laminates. 
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TABLE 1 


MATRIX PROPERTIES VS TEMPERATURE 
(061*0 ALUMINUM ALLOY 


TEMPERATURE, 

1 

MODULUS, 
10® PSI 

POISSON'S 

RATIO 

THER COEFF EXP, 
10** 

YIELD 

STRESS, 

PSI 

0 

113 

a 310 

12.1 

8500 

100 

ILO 

.330 


8500 

200 

107 

.350 


8500 

300 

ia4 

.345 


8000 

400 

lai 

.335 

14.7 

7000 

500 

9.9 

.330 

15.2 

4000 

600 

9.8 

.330 

15.8 

3000 

700 

9.3 

.330 

1&5 

2000 

800 

&8 

.345 

17.2 

1000 

900 

ao 

.350 

17.9 

0 


7.0 

.350 

ia6 

0 


TABLE 2 

MATRIX STRESS/STRAIN PROPERTIES 
0061*0 ALUMINUM ALLOY 


STRAIN, 

% 

— 

STRESS, 

PSI 

MODULUS, 

PSI 

POISSON'S RATIO 

0 

0 

ILOx# 

a366 

.024 

2 590 

ia5 

.336 

.030 

5 260 

9.8 

.334 

.078 

7920 

&2 

.332 

.098 

9 250 

3.6 ^ 

.392 

.44 

10 500 

l7xW 

.409 

12 

11900 

a&(i(r 

.400 

L5 

11900 

9.9xl0j 

.397 

5.0 

11900 

LOxlO* 

.450 















TABLE 3 


i 

i 

I 

j 


I 

I 

I 

I 
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MATRIX SHEAR STRESS-STRAIN 


PROPERTIES 


6061-0 ALUMINUM AaOY 


STRAIN, 

% 


MODULUS, 
106 PSI 

0 

0 

3.27 

.04 

1150 

3.44 

.07 

2450 

3.57 

.11 

3760 

3.56 

.15 


3.32 

.19 

6430 

2.03 

.33 

7810 

.867 

.50 


.100 

100 

8300 

.010 


TABLE 4 

CUMULATIVE FIBER FRACTURE 


STRESS LEVEL 
(KSI) 

FIBERS FRACTURED 
(PERCENT) 

200 

0 

250 

3 

300 

5 

350 

6 

400 

15 

450 

30 

475 

50 

500 

72 

525 

90 

530 

95 

550 

100 











Figure L - Flow chart ol nonlinear laminate analysis for matal matrix filiar composites. 
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Figurt 2. * Cwnptrison 01 mHturtd and predicted results. 
Boron/aluminum unidirectional composites. (4 mil diam. 
fitwr. 6061-AI, 0, SO fiber volume ratio!. 



Figure 1 - Nonlinear intraUmnar sfiear stress-strain curve of boron/ 
aluminum. 



Fifurt 4. - Comparison o( mtasurtd and prtdicltd rtsults. Boron/alumi- 
num ( UtlJi, laminato loaded -B0° to 0^*ply. (4mildiam. fiber, iObl- 
Al, 01 Sofiber volume ratio). 
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Figure S. - OP-ply ftress-stnln behavior. BororValumlnum I MilSlj 
laminate. (40inlldlam. ftter. M61*AI, dS fiber volume ratlol. 







Figure 6. - IS^-ply stress-strain behavior. Boron/aluminum [ UtlSlj 
laminate. (40mUdiam. fiber, 6061-AI, OL SO fiber volume ratio). 



Figure 7. - oP-jHy cumulative stresses In boron/ 
aluminum { (M±15|, laminate. (4. 0 mil diam. 
fiber. 6061 -aI, (X5 fiber volume ratio). 
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Figure & - Nonlinear laminate analysis results compared with measured data. 



Figure 9. - OP-ply stress-strain behavior. Boron/aluminum [0^±15jj 
laminate. (40mildiam. fiber, 6061-AI, OLM fiber volume ratio). 










